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Abstract 

The photo-disintegration of cosmic ray nuclei by starlight leads to the production of sec¬ 
ondary antineutrinos. We have assumed that the flux of the ultrahigh energy cosmic ray nuclei 
near the Galactic plane region is the same as that observed near the earth and calculated the 
antineutrino flux produced from their photo-disintegration. The IceCube detector has mea¬ 
sured the neutrino/antineutrino flux in the TeV-PeV energy range. Our calculated secondary 
antineutrino flux in the energy range of 10-100 TeV is found to be much less compared to the 
flux detected by the IceCube collaboration. The upper limit on the intensity of the radiation 
field in the extragalactic medium is much lower than that near the Galactic center. If we extend 
our formalism to the extragalactic medium the contribution from the photo-disintegration of 
ultrahigh energy cosmic ray heavy nuclei remains insignificant due to their very low flux. 

Keywords: high energy cosmic rays, photo disintegration 

1 Introduction 

The cosmic ray detectors of CASA-MIA [1], Tunka [2], GAMMA [3] and KASGADE-Grande [4] 
have measured the cosmic ray flux from the knee region at around 10^^ eV which mostly contains 
iron nuclei. This is known as the “iron knee”. Beyond the knee there is a change in the composition 
of the cosmic rays from heavy to light nuclei. 

The radiation held produced from the stars is most intense near the Galactic center region where 
the density of the stars is the highest and gradually decreases with the distance [5] from this region. 
The cosmic rays emitted from the Galactic and extragalactic sources interact with the interstel¬ 
lar radiation and matter during their propagation through the interstellar medium. Secondary 
neutrinos and gamma rays are produced from their interactions [6, 7, 8, 9]. The pure hadronic 
interactions have negligible contribution to the observed neutrino events in the IceGube detector 
[10, 11, 12, 13]. 

The photo-disintegration of cosmic ray heavy nuclei has been studied earlier analytically [15, 16, 17] 
and with Monte Garlo simulations [18]. More recently it has been suggested that the antineutri¬ 
nos produced in the photo-disintegration [19] of some very high injected flux of ultrahigh energy 
cosmic rays (UHEGRs) may explain the flux observed by the IceGube neutrino detector. In photo¬ 
disintegration protons and neutrons are produced almost equally. The neutrons decay to protons 
and antineutrinos. In this paper it has been assumed that 1 to 10% of the high energy protons 
produced in the decay of the high energy neutrons reach us due to magnetic shielding and their 
flux is the same as the observed cosmic ray proton flux in the energy range of 10®’® to 10®’® GeV. 
Parametrizations of the all particle cosmic ray spectrum are given in [20] for the following three 
population models (i) Hillas model (ii) global fit model. We have used their parametrizations of 
the diffuse cosmic ray flux to calculate the antineutrino flux produced in the photo-disintegration 
of cosmic ray heavy nuclei. According to the Hillas model the Galactic cosmic ray spectrum ends 
at the knee which mostly originates from supernova remnants. The contribution to the diffuse 
cosmic ray flux from the extragalactic sources is significant at the ankle. This scenario is based on 
the amplification of the magnetic helds in non-linear diffusive shock accelerations which determines 
the maximum energy of the cosmic rays produced in Galactic SNRs. In this model at least three 
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populations of particles are needed to explain the observed cosmic ray flux. Many more populations 
may be introduced to explain the cosmic ray spectrum in much more detail. In the global fit model 
the fluxes of the individual cosmic ray nuclei measured by the CREAM experiment [21] are well 
explained. This model is also consistent with the “iron knee” observed by the KASCADE-Grande 
at 10^^ eV [22]. In each population all the particles have the same maximum rigidity. Another 
important aspect is the higher energy populations can significantly contribute to the cosmic ray 
flux at lower energy. For more details on these models the readers may see [23]. 

Although, the radiation field depends on the Galactocentric radius R and the altitude above the 
Galactic plane z, a uniform background radiation field is assumed in this work near the Galactic 
plane region. The strength and spectral distribution of this radiation field are assumed as given 
for R=0,4kpc and z=0 in Fig.l.(bottom panel) of [5]. We have denoted the two cases by IRl and 
IR2. In the first case (IRl) it is assumed that the Galactic plane region has a uniform radiation 
field, the same as the one given for R=0,z=0 in their paper. IR2 is the case corresponding to the 
radiation field equal to that given for R=4kpc,z=0. 


2 The Cosmic Ray Flux Near the Galactic Plane: 


The cosmic ray flux and its composition near the Galactic plane region are unknown. Even near 
the earth the composition of the UHEGRs is not yet known. The measured flux near the earth has 
been fitted with the Hillas model and the global fit model earlier [20]. The all particle spectrum 
(GeV“^ m~^ sec~^ sr“^) in the three population models is 


dNA^EAj 
dEAi dS dt dQ 


= ^ X exp 


EAi 1 


( 1 ) 


We have used the above cosmic ray flux in this work. The subscript i = 1, 5 runs over the standard 
five groups of particles p. He, GNO, Mg-Si and Fe. The three populations are denotd by j = 1,3. 
Area has been denoted by S and charge of a nucleus by Zi. For GNO and Mg-Si we have used the 
mean values of their charges. The values of the parameters are given in Table 1 and Table 2 for 
the Hillas model and the global fit model respectively. Note that in eqn.(l) the energy EAi is the 
energy of the nucleus for the group, EAi ~ En where Aj is the average mass number for the 
j^th gj-Q^p ig energy of each nucleon. In the Hillas model the fluxes of the heavy nuclei 

are higher compared to the global fit model. The population 3 in the global fit model contains only 
protons and iron nuclei whereas in the Hillas model proton. He, GNO, Mg-Si and iron fluxes are 
almost equal. In the next section we do not use the subscript i anymore to write the cosmic ray 
flux for the individual groups. 


3 Secondary Antineutrinos from the Galactic Plane Region: 


The cross-section of photo-disintegration [24] for a medium or heavy nucleus of mass number A in 
the rest frame of the nucleus is 


— a^A 




(e*^ — Cg )2 -|- (e*A)2 


(2) 


where the photon energy in the nuclear rest frame e* is below 30 MeV. Above 30 MeV the cross- 
section is energy independent with a smaller value of A/8 mb. The values of the constants are 
as follows, cross-section cjoa = 1.45 A mb, the central value of GDR (Giant Dipole Resonance) 
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eQ = 42.65 yl MeV for mass number of nuclei yl > 4 and width of the GDR A = 8 MeV. The 
rate of photo-disintegration {Rph) is calculated for the background radiation density per unit energy 
(IR) The Lorentz factor of each nucleon of mass rUn and energy En is = Ea/{A mn)- 


Rph — 






^th 


e*aA{e*)de* 


‘/An 


dn{x) dx 
dx 


(3) 


The threshold energy of photo-disintegration in the rest frame of the nucleus is €th ~ 2 MeV and 
the maximum energy of the photons in the background radiation held is €max- 
Eq.(3) is simplihed after including the expression for the total cross-section of photo-disintegration. 
The cross-section corresponding to the GDR can be approximated by a delta function [17]. 

aA{e*) = naoAY^ie* - e'o) (4) 


The expression for the rate of photo-disintegration in the case of GDR is simplihed to 




■ph,GDR 


CTraoA e'o A 


dn(x) dx 


47n A[,/ 27 n dx x^' 


The constant cross-section of photo-disintegration above 30 MeV gives 


R 


Ac r'^'yn^r} 


ph,const 


IGjn hoMeV 


e*de* 


j-emax dn{x) dx 
/(e‘/ 27 n) dx x 2 


(5) 

( 6 ) 


The contribution from the constant cross-section above 30 MeV is found to be low. We have 
considered a disc shaped region centered at the Galactic center of radius R = lOkpc and height 
z = 0.5 kpc having uniform intensity of radiation. The rate of photo-disintegration is highest for 
Fe and decreases with decreasing A as shown in our Fig.l. The spectrum of neutrons [GeV“^cm“^ 
sec“^] produced in photo-disintegration is 


dNn^ {Ens ) 

dEns dt dV 


0.5Rph{EnJ 


dNAjE^J 
dEns dV 


(7) 


Both neutrons and protons are stripped out from the nuclei. We have assumed 50% of the stripped 
nucleons are neutrons and denote them by n^. The energy of the stripped neutrons is assumed to 
be the same as that of the parent nucleons. The steady state density of cosmic ray nuclei is 
expressed in per nucleon energy in eqn.(7). This density flux (GeV“^ cm“^) is 


dNA{En) ^ 47r AdNAjEA) 

dEn dV c{cmfsec) dEA dS dt dQ 


The neutron flux [GeV ^ cm ^ s ^ sr is calculated for the disc shaped region. The geometrical 
correction is done following the formalism discussed in paper [8]. 


Jus (d^ris ) 


^ r dV dNnsjEns) 
Dg J 4 tt dEns dt dV 


(9) 


r is the distance of the site of neutron production from the earth. The geometry of our Galaxy is 
symmetrical around the Galactic center. We have calculated the effective radius of the halo with 
respect to the observer on earth assuming a homogeneous spherical halo centered at the earth. This 
effective distance or radius is obtained after integrating over r, Re// = / dV/{4'K r‘^) =1.7 kpc as 
discussed in [8]. 
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Figure 1: The positions of the Super Nova Remnants (SNR), Low and High Mass X-ray Bina¬ 
ries (LMXB, HMXB) from http://vizier.cfa.harvard.edu, Central Molecular Zone (CMZ), Galactic 
Molecular clouds (M-cloud) from https://www.cfa.harvard.edu are plotted in this skymap. The 
IceCube detected events (diamond shaped points) [13] and the Fermi Bubbles (two patches) are 
also shown [25]. 


The positions of the Galactic astrophysical objects are shown in the skymap in Fig.l. The solid 
lines show the region \b\ < 10° which includes most of the photon emitting objects. We have 
selected an even smaller region \b\ < 5° to define the Galactic plane region in the present work. In 
this region the intensity of starlight is expected to be very high. The solid angle subtended by this 
region to us is =1 sr. 

The neutrons decay to antineutrinos with average energy = 0.48 MeV in the rest frame of the 
neutrons. The secondary neutrons have a flux given by eqn.(9). The antineutrino flux (GeV“^ 
cm“^ sec“^ sr“^) produced from the decay of neutrons is 


ME^) 


rUn 

2eL 


clEn 


Eus 


'Jfis i,En 


( 10 ) 


as given in eqn.(25) of [17] where E^^^ is the maximum energy of the neutrons decaying to an- 
tineutrinos. The antineutrinos may take energy between 0 and We have assumed that the 

antineutrinos on average have energy to simplify the calculations. The number of neutrons 

is equal to the number of antineutrinos if all the neutrons decay before reaching us. We apply the 
conservation of the total number of neutrons and antineutrinos to get the antineutrino flux. The 
antineutrino flux (GeV“^ cm~^ sec“^ sr“^) produced from the decay of neutrons of energy En^ and 
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Parameters 

P 

He 

GNO 

Mg-Si 

Ee 

Pop 1: 

Uij 

7860 

3550 

2200 

1430 

2120 

Rc = 4PV 

lij 

1.66 

1.58 

1.63 

1.67 

1.63 

Pop 2: 

ttij 

20 

20 

13.4 

13.4 

13.4 

Rc = 30PV 

lij 

1.4 

1.4 

1.4 

1.4 

1.4 

Pop 3 

Uij 

1.7 

1.7 

1.14 

1.14 

1.14 

Rc = ‘lRN 

lij 

1.4 

1.4 

1.4 

1.4 

1.4 


TABLE 1 : Hillas Model: normalisation constants aij and spectral indices 7 jj and cut-offs in rigidity 
Rc (Gaisser, Stanev& Tilav 2013). 



Parameters 

p 

He 

G 

O 

Fe 

Pop 1: 

aij 

7000 

3200 

100 

130 

60 

Rc = 120TV 

lij 

1.66 

1.58 

1.4 

1.4 

1.3 

Pop 2: 

aij 

150 

65 

6 

7 

2.3 

Rc = 4PV 

lij 

1.4 

1.3 

1.3 

1.3 

1.2 

Pop 3 

aij 

14 




0.025 

Rc = 1.3 EV 

lij 

1.4 




1.2 


TABLE 2 : Global fit Model: normalisation constants aij and spectral indices 7 ij and cut-offs in rigidity (Gaisser, 
Stanev & Tilav 2013). 


Lorentz factor is 

J,(E,) = ^JnAEnJ ( 11 ) 

where Ep = '^ns^'n- We note that the antineutrino fluxes calculated using eqn.(lO) and eqn.(ll) 
differ by a factor of 2““^/a where a is the spectral index of the stripped neutron spectrum. For 
a in the range of 1 to 4 one gets accurate results within a factor of 2 using our eqn.(ll). Our 
calculated antineutrino flux is found to be very low as shown in Fig.3. 

4 Photo-Disintegration of UHECRs in the Extragalactic Medium 

The formalism discussed in the earlier section is extended to calculate the antineutrino flux produced 
in the photo-disintegration of ultrahigh energy cosmic ray nuclei during their propagation from 
extragalactic sources to the earth. 

It has been suggested recently that the photo-disintegration of some very high injected flux of 
extragalactic cosmic ray nuclei of energy between 10®'^ GeV/nucleon to 10®'^ GeV/nucleon by the 
far infrared background at 10 meV may explain the neutrino/antineutrino flux detected by the 
IceGube detector in the energy range of 10®'® to 10®'^ GeV [19]. The flux of the injected or parent 
ultrahigh energy cosmic ray nuclei is unknown in this paper. When the ultrahigh energy cosmic 
ray nuclei are photodisintegrated high energy protons and neutrons are produced. The high energy 
neutrons decay to high energy protons and antineutrinos. The author has assumed that only 1 to 
10% of the high energy protons reach us due to shielding by the extragalactic magnetic field and 
their flux is the same as the observed proton flux in the energy range of 10®'^ to 10®'® GeV. The 
antineutrino flux from the decay of the high energy neutrons is calculated in this paper by assuming 
that this neutron flux is 10 to 100 times higher than the observed proton flux in the energy range 
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of 10®'® to 10®'® GeV. Due to this assumption an antineutrino flux of ~ 10“® GeV cm“^ sec“^ sr“^ 
is estimated in [19]. 

The approach follo'wed in our paper is completely different. We have used the steady state 
or observed flux of ultrahigh energy cosmic ray nuclei 'with energy 10®'® GeV/nucleon to 10®'® 
GeV/nucleon to calculate the flux of the high energy neutrons produced in photo-disintegration. 
The antineutrino flux produced in the decay of these neutrons is subsequently calculated. Thus the 
normalisation used in our paper does not depend on any assumption or magnetic shielding. 

We have assumed a spherical cosmological volume within which the cosmic ray sources are dis¬ 
tributed isotropically. The isotropic UHECR flux within this volume is assumed to be similar to 
that predicted in the Hillas model. The heavy nuclei are photo-disintegrated in the ambient radia¬ 
tion field and antineutrinos are produced through the decay of the UHECR neutrons. At 10 meV 
the intensity of the radiation field is 1 eV/cm® near the Galactic center according to the radiation 
fields (IRl and IR2) calculated in [5]. The intensity of the extragalactic background light at this 
energy has un upper limit of 0.005 eV/cm® from observations on gamma ray sources [26], which is 
200 times lower than that at the Galactic center. The energy of the ambient low energy photons is 
boosted in the nuclear rest frame by the Lorentz factor of the nucleus which disintegrate the cosmic 
ray nucleus. The contribution from the GDR cross-section is maximal for incident photon energy 
in the energy range of 10-30 MeV in the nuclear rest frame. The CMB is present throughout the 
interstellar medium. However, for photo-disintegration by GMB the UHECRs have to be of energy 
more than 10^® GeV/nucleon and the secondary antineutrinos would have energy more than 10® 
GeV. At present we do not have any observational result on the neutrino/antineutrino flux at this 
energy. 

Our Eqn.(lO) has been modified to find the flux in the observer’s frame 

Tob/pob'i _ J_ f dNn^jEn^) . . 

dEnJtdVc ' ^ ^ 


The comoving distance and differential volume are dc and dVc respectively. The effective radius 
of the cosmological volume centered at the observer on earth is Reff = / dVc /Within 
this volume the flux of the ultrahigh energy cosmic ray nuclei is assumed to be the same as their 
observed flux near the earth. The densities of the radiation field and the ultrahigh energy cosmic 
rays are assumed to be uniform within this volume. If the redshift corresponding to the distance 
Reff is Zeff then the energy and time in the observer’s frame and at redshift Zg// are related as 
E°^ = E/{1 + Zeff), t°^ = t{l + Zeff)- The antineutrino flux per unit energy unit time in the 
observer’s frame is 


job f zT'ob \ 




(13) 


Eqn.(12) and Eqn.(13) are used to estimate Reff- The AGN which are bright in X-rays and 
potential sources of extragalactic cosmic rays, are located near redshift 2 to 3 [27]. Assuming 
Zeff ~ 2 the correction to the energy of the antineutrinos in the observer’s frame can be neglected. 
To produce the observed flux of 6.7 x 10“®® GeV“®cm“^ sec“® sr“® near 10® GeV [28] the required 
effective radius is found to be unrealistically high ~ 50000 Gpc. As the density of the UHEGR 
nuclei is very low a large effective radius or volume is required to explain the flux detected by the 
IceGube detector. 


5 Discussion: 

We have used the three population models (i) Hillas model (ii) global fit model of cosmic rays to 
calculate the flux of antineutrinos produced in the photo-disintegration of cosmic ray heavy nuclei 
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Rate Plot 



Figure 2: Rate of Photo-disintegration of cosmic ray iron and oxygen nuclei plotted for the assumed 
radiation fields IRl (R=0,z=0) and IR2 (R=4kpc,z=0). 

near the Galactic plane region. The CNO, Mg-Si and Fe cosmic ray nuclei are included in our study. 
The radiation field is assumed to be uniform within the disc shaped region of radius 10 kpc and 
height 0.5 kpc centered at the Galactic center. The spectral intensities of the field corresponding 
to R=0,z=0 and R=4kpc, z=0 given in Fig .1. (bottom panel) of the paper by [5] have been used 
in the present work. 

The cross-section of photo-disintegration has a resonance below 30 MeV and above 30 MeV it is a 
constant with a smaller value. We have calculated the rates of photo-disintegration of oxygen nuclei 
in the radiation field IRl for the resonance and constant values of the cross-section. The contribu¬ 
tion from the constant cross-section is found to be low. Our calculated rates of photo-disintegration 
for oxygen and iron nuclei (Rph) are plotted in Fig.2. for the radiation fields defined at R=0,z=0 
(IRl) and at R=4kpc,z=0 (IR2). The energy dependence of the rate of photo-disintegration can be 
explained from the energy dependence of the background radiation field. The higher energy cosmic 
ray nuclei can be photo-disintegrated by the lower energy photons in the radiation field. 

If we increase the halo size the effective radius increases (please see Table 1: of [8]). The antineutrino 
flux is directly proportional to the effective radius so it also increases. 

The IceGube collaboration has measured the isotropic flux of neutrinos and antineutrinos in the 
TeV-PeV energy range. In the energy range of 25 TeV to 2.8 PeV their measured flux is = 

6.7 X 10“®(Fl,y/100TeR)“°'® GeV cm“^ sec“^ sr“^ [28]. Our calculated fluxes from the Galactic 
plane region shown in Fig.3. are many orders of magnitude lower than the flux detected by the 
IceGube detector. The average length traversed by the neutrons before decaying to antineutrinos 
increases with energy. At 10^ GeV this length is 1 kpc. We do not expect antineutrinos from the 
Galactic plane region at very high energy as their parent ultrahigh energy neutrons will not decay 
before reaching us. 

After showing that the antineutrino flux produced in the Galactic plane region is very low we 
have extended our formalism to the extragalactic medium. Assuming the sources are located at a 
effective redshift 3 we have calculated the antineutrino flux from photo-disintegration of the flux 
of UHECR nuclei in the global fit model by the extragalactic background light. The upper limit 
on the intensity of extragalactic background light has been used from [26] in our calculation. Our 
calculated antineutrino flux is shown in Fig.4. In this case also the contribution to the observed 
flux from the photo-disintegration of UHECR nuclei is very low. 
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Flux of Antineutrinos 



Figure 3: Antineutrino fluxes from the Galactic plane region calculated in the energy range of 10 
TeV to 100 TeV for the Hillas model and the global fit model assuming uniform radiation fields for 
the two cases IRl and IR2. 


Flux of Extragalactic Anlineutrinos 



Figure 4: Maximum flux of antineutrinos calculated for effective redshift of sources Zeff = 3 and 
maximum extragalactic background light 0.005 eV/cm^ [26]. The global fit model of the UHECR 
flux has been used. 

6 Conclusion: 

The three population models are used to calculate the antineutrino flux produced in the photo¬ 
disintegration of CNO, Mg-Si and Fe nuclei near the Galactic plane region. In all cases it is found 
that the antineutrino flux is much less than the flux measured by the IceCube detector. The 
antineutrino flux is directly proportional to the effective radius of the spherical region within which 
the cosmic ray nuclei are photo-disintegrated. The formalism discussed for the Galactic plane 
region has been extended to study the contribution from the photo-disintegration of cosmic ray 
nuclei outside our Galaxy. In this case also the antineutrino flux is very low compared to the 
observed flux due to the low flux of the parent UHEGR nuclei. 
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